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ABSTRACT 

This paper Investigates a design method of 
an active control bearing system with only veloc- 
ity feedback. The study provides a new quasl- 
modal control method for a control system design 
of an active control bearing system In which 
feedback coefficients are determined on the basis 
of a modal analysis. Although the number of sen- 
sors and actuators Is small, this quasl-modal 
control method produces a control effect close to 
an Ideal modal control. 

1. INTRODUCTION 

There are two approaches to reduce the vibra- 
tion of rotating shafts. One Is the approach 
using passive elements and the other Is the case 
using active elements. 1-13 The studies of the 
former based on damped flexible supports have 
been done by many authors. However, It Is diffi- 
cult to obtain desirable optimum tuned conditions 
In actual rotating machinery. Conversely, It Is 
very easy for an active vibration control method 
to obtain desirable optimum values because the 
support conditions can be varied by only feedback 
coefficients. The first author proposed an active 
control bearing wherein the bearing housings are 
actively controlled by using a state feedback 
control method. 12 * 13 The experiment for a two 
degree of freedom system proved that active con- 
trol bearings are effective for vibration control 
of rotating shafts. However, for a multidegree 
of freedom system. It Is necessary for a control 
system with state feedback to get complete state 
variables of displacements and velocities for a 
rotor system. A control system Is thus compli- 
cated. Therefore, In the case of a multidegree 
of freedom system such as a multi-bearing and 
multi-disk system, active control bearings with 
only velocity feedback are to be desired. This 
paper Investigates a design method of an active 
control bearing system with only velocity feed- 
back. The study provides a new quasl-modal con- 
trol method for a control system design of an 
active control bearing system In which feedback 
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coefficients are determined on the basis of a 
modal analysis. The features of the new quasl- 
modal control method proposed In this paper dif- 
fer depending on the number of measurements and 
the number of controlled modes. There are three 
cases. The first case Is when the number of 
measurements of vibration velocity Is equivalent 
to the number of controlled modes. The second 
case when the number of measurements Is less than 
the number of controlled modes. The third case 
Is when the number of measurements Is more than 
the number of controlled modes. The method does 
not depend on the number of active control 
bearings. 

On the basis that It Is sufficient for 
active control of a system to provide critical 
damping on each mode, the control system Is 
assembled independently for each mode. Simula- 
tions of the active control of a three bearing 
and three disk rotor system are carried out to 
verify the efficiency of this method. Although 
the number of sensors and actuators Is small, 
this quasl-modal control method produces a con- 
trol effect close to an Ideal modal control or an 
optimal state feedback control. 

2. QUASI-MODAL VIBRATION CONTROL METHOD 

2.1 Ideal Modal Control System Design Method 

First, a control system design method Is 
described In general. The following equation Is 
considered. 

Hnin * £ n X n * K n X n - P(t) ♦ U(t) (1) 

where M n , C n , and K n are mass, damping, and 
stiffness matrices. X n Is a general nodal dis- 
placement vector. P(t) Is a general nodal unbal- 
ance force vector and U(t) Is a general nodal 
control Input vector. Equation (1) Is written with 
n degree of freedom system. For X n of Eq. (1), 
the following transformation Is performed. 

X n = Ta or a = T”^X n (2) 

where T Is a modal matrix assumed to be normal- 
ized. With Eq. (2), Eq. (1) Is transformed as 
follows : 

I T 5nlK f I T £nli ♦ I T K n Ia = T T P(t) + T T U(t) (3) 
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where T T M n T Is a unit matrix, T T K n T Is a 
diagonal frequency matrix, and 


11 


21 


I T £nI - 


C 12 * * * C 1 n 

C 22 * * * C 2n 


nl L n2 ' 


. C 


If constant damping coefficients (proportional 
damping) are assumed. 


where 


G = (T T ) - ^ FT - ^ 


( 11 ) 


In n degree of freedom system, If the number of 
control Inputs or number of velocity measurements 
Is less than n, this control system Is not a 
true modal control system. This paper names such 
a control system a quasl-modal control system. 

For example, for two control force Inputs and 
three vibration measurements, the quasl-modal 
control system Is expressed as follows: 


U * G'Xn 


( 12 ) 
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Therefore Eq. (3) Is reduced to the following 
equation. 

a + ca +■ fi 2 a = T T P(t) * T T U(t) (6) 


where ST symbolizes a frequency matrix, and c 
Is a principal damping matrix. Thus Eq. (1) Is 
reduced to Eq. (6) separated on each mode and 
uncoupled. It Is very easy to determine the 
optimal control Inputs to each mode In the above 
expressed modal domain. The Input T T U(t) Is 
taken to be a portion of critical damping times 
the modal velocity a. That Is, 


F = 


- 2f n 


T 

T U = Fa 
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(7) 


( 8 ) 
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and Ci, C 2* Cn are modal damping ratios, 

or percentage of critical damping, and w-j , <*> 2 , 
..., w n are undamped critical speeds. 

If modal damping ratios Inherent In the sys- 
tem (obtained by measurement) are very small, 

let Ci = 0.7 (1 = 1, 2 n) In Eq. (8). If 

the system damping Is not negligible, then the 
values of Cl In Eq. (8) should be taken as 0.7 
minus the measured value. The control Input, 
from Eq. (7) , Is then 


T -1 • 
U = (T* ) 'Fa 


( 9 ) 


or, using the definition of a (Eq. 2), 

U = GX n (10) 
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where G* Is a modified feedback coefficient 
matrix."* Determining G' Is the subject of the 
next section. 

2.2 Modification of Feedback Coefficient for 
Quasl-Modal Control 

Equation (11) provides feedback coefficients 
for the case when vibration velocities are meas- 
ured at every station. If this Is not the case, 
satisfactory control Is unlikely using Eq. (11). 
Therefore for a quasl-modal control system, It Is 
necessary to modify the value of G depending on 
the number of measurements and positions of meas- 
urements (probe locations) In order to produce 
the maximum of control efficiency. The assump- 
tions are as follows; (1) the number of active 
control bearings and their positions are given, 

(2) the vibration from the first mode to the s th 

mode (s < n) Is to be controlled, and (3) the 
positions of measurements (probe locations) are 
unconstrained. The modification of G differs 
depending on the number of measurements and the 
number of controlled modes. They are divided 
Into three cases depending on the relationship of 
the number of measurements (k) to the number of 
controlled modes(s): (1) k = s, (2) k < s, and 

(3) k > s. 

2.2.1 Case when the number of measurements 
of vibration velocity Is equivalent to the number 
of controlled modes , (k = s). As the first step, 
write the G matrix with zero rows where there 
are no active control bearings. 

locations of measurements 
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From Eq. (14), the locations of active control 
bearings are positions In which row vectors in 
the matrix G are not zero. 

For the first mode, the control force 
coefficients supplied to each active control 
bearing may be computed by using G from 
Eq. (14). 
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where [ty|, t2i » ■ • • , t n ] is the eigen vec- 
tor for the first mode In the modal matrix T. 
Namely, this vector multiplied by an angular 
velocity co and modal amplitude Is the first 
modal component of the vibration velocity vec- 
tor. On the right side of Eq. (15), v a -| , v^ , 
...» v e i are the control force coefficients sup- 
plied to positions a, b, . . ., and e to 
actively control the first mode vibration. In 
the same manner for the other modes, the follow- 
ing relations are obtained. 
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Combining these relations, 


The right side of Eq. (17) shows the control force 
coefficients supplied on each mode to each posi- 
tion of active control bearings. The control 
forces applied to any one bearing are Independent 
of the forces applied to the other bearings. With 
respect to the active control bearing at position 
a, It Is the most desirable to supply control 
forces of v a -j for the first mode, v a 2 for the 
second mode, . . ., and v as for the s th mode. 
Therefore, even If the number of vibration veloc- 
ities measured are less than n, we only have to 
modify the elements In the matrix G to satisfy 
the right side of Eq. (17). 

As the number of measurements (k) equals the 
number of controlled modes (s) from assumption, 
replacing n with s and writing Eq. (17) again 
after omitting zero elements, we obtain the fol- 
lowing expression. 
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(18) 


or 

G'T = V 


(19) 


where the dimensions of the matrices G' and V 
are r x s and the matrix T has the dimension 
s x s. From Eq. (19), we have 
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or 


G' = VT- 1 (20) 

The matrix G' Is the modified feedback coeffl 
clent matrix to be used Instead of the matrix 
G. In Eq. (11) G' Is the one and only solution. 

The number of active control bearings Is 
equivalent to the number of row vectors of the 
matrix G' In Eqs. (20) and (21), namely r. 
Whether r > s or r < s, the matrix G' Is cer- 
tainly determined. However, the eigen value 
assignments and the unbalance responses actually 
depend on the number of active control bearings 
as shown later. In this section, an optimization 
of control Inputs Is carried out concerning how 
to supply the control forces to active control 
bearings assigned. Even If the number of measure- 
ments of vibration velocities Is at most equiva- 
lent to the number of controlled modes, It Is 
possible to achieve the control effect similar to 
the case when the number of measurements Is n. 

It Is important to note that this control effect 
does not depend on the location of measurements. 
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2.2.2 Case when the number of measurements 
of vibration velocity Is less than the number of 
controlled modes , (k < s). In this case, Eq. (18) 
appears as 
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where the matrix G* has the dimension r x k 
(k < s), the matrix T has k x s and the matrix 
U has r x s. Thus T becomes a rectangular 
matrix and cannot be Inverted. Therefore, the 
number of controlled modes Is reduced to k to 
be able to get the Inversion of the matrix. 

After this reduction, matrices T and G" are 
obtained as follows: 
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s the following 

consideration 


Accordingly, 


for the modes eliminated from a control object. 


Assuming that the modes eliminated are 1 th and 
j th, we estimate errors as regards the 1 th and 
the j th modes between the control force deter- 
mined by Eq. (25) and the optimal control force 
In section 2.2.1 using the matrix G' (Eq. (18)) 
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where v a1 , v b i v e1 and v a j, v b j, - - 

v e j of the right side are control forces supplied 
for the 1 th and the j th modes eliminated. 
Accordingly, k positions of measurements have to 
be optimized so as to satisfy the following per- 
formance Index. 
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The dynamic characteristics sensitivity depends 
on k positions of measurements In this case. 

In general. It Is Impossible to realize J = 0; 
however It Is possible to make J small. It Is 
very Important to find out the optimal positions 
of measurements. Moreover, It has to be also 
evaluated which mode can minimize the performance 
Index J. 


2.2.3 Case when the number of measurements 
of vibration velocity Is more than the number of 
controlled modes , (k > s). This case is k > s 
In contrast to section 2.2.2. We only have to 
determine the feedback coefficients by Increasing 
controlled modes as to k = s. 

The three cases above mentioned are summa- 
rized as follows: (1) In the first case of 

k a s, unique feedback coefficients are deter- 
mined Independent of measuring positions and the 
number of active control bearings, (2) In the 
second case of k < s, the feedback coefficients 
are not determined such as the first case. 
Accordingly, the number of controlled modes must 
be decreased so the number of controlled modes 
agrees with the number of measurements. Then the 
feedback coefficients can be determined. After 
this, on the basis of a performance Index for 
uncontrolled modes, the positions of measurements 
must be selected so as to minimize the perform- 
ance Index. These positions of measurements are 
the best positions to measure and the feedback 
coefficients In this case are best for vibration 
control. (3) In the third case of k > s, as In 
the second case, feedback coefficients are not 
determined. In this case, contrary to the second 
case, the number of controlled modes must be 
Increased . 

The control effect In the case where veloc- 
ities are measured at all positions can be real- 
ized by measuring velocities at only s positions 
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minimize the performance Index. It Is possible 
to control a vibration of rotating shaft from the 
first mode to the third mode as shown In Fig. 5 
(a) and (b). The first mode Is evaluated for the 
performance Index J In this example. 

4. CONCLUSION 

According to the results of the simulations, 
If the optimal feedback coefficients are chosen, 
the unbalance vibration up to the third mode can 
be sufficiently controlled by means of only two 
active control bearings and measurements of vibra- 
tion velocities at only two locations. The unbal- 
ance amplitude can be reduced to less than the 
center of gravity eccentricity. In this case, It 
Is very Important to choose the positions of meas- 
urements to avoid Instability. 

The quasl-modal control method for active 
control of rotor vibrations proposed In this paper 
results In near maximum control efficiency using 
the minimum number of active control bearings and 
the minimum measurements of vibration velocities. 
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In the first case and the third case. In all 
three cases, the greater the number of active 
control bearings, the better the response of 
vibration control will be. 

3. MODEL OF ROTOR-BEARING SYSTEM AND SIMULATION 

3.1 Rotor Model and Equation of Motion 

Simulations of a three bearing and three 
disk rotor system are carried out In this chap- 
ter. The rotor model and Its undamped natural 
modes are shown In Fig. 2. The equation of 
motion Is derived by a finite element method. In 
this case, the equation of motion Is a fourteenth 
order matrix differential equation because this 
rotor model Is divided Into six elements. Now, 

It Is assumed that there are three active control 
bearings at positions 1, 3, and 6 and this 
system is uncoupled between x and y directions. 
The vibrations from the first mode to the third 
mode are to be controlled and It Is possible to 
measure vibration velocities at any position. 

For simplicity, only measurements In the x direc- 
tion are considered. 

3.2 Case When the Number of Measurements of 
Vibration Velocity Is Equivalent to the 
Number of Controlled Modes 

This section describes the case where the 
number of measurements of vibration velocities Is 
three. In order to Investigate the stability of 
the rotor-bearing system, a complex eigenvalue 
analysis Is carried out and the eigenvalues are 
shown on a complex plane as an eigenvalue assign- 
ment. The unbalance responses In these cases are 
also shown. Figure 3 (a) and (b) shows the eigen 
value assignment and their unbalance responses. 

In the case of three Inputs and three outputs 
like this one, It Is observed that the unbalance 
responses are almost the same independent of the 
positions of measurements. The unbalance 
responses are similar to the case In which all 
velocities at all positions are measured. This 
section also shows the responses when the number 
of active control bearings Is changed. These 
results prove the summary In chapter 2. 

3.3 Case When the Number of Measurements of 
Vibration Velocity Is Less than the 
Number of Controlled Modes 

This section Illustrates the case of three 
control Inputs and two measurements as the second 
case In chapter 2. The control characteristics 
depend on the mode of the performance Index. The 
positions where the performance Index Is minimum 
should be selected for measuring positions of 
vibration velocities. These are shown In Fig. 4 
(a) and (b) as typical example. These show the 
case In which the performance Index J Is based on 
the third mode. These figures Indicate that 
measuring positions x-| and X 2 are the best 
positions In this case. Since some cases become 
unstable In such a system with three Inputs and 
two outputs, care Is required In the selection of 
measuring positions to avoid having unstable 
modes . 

Lastly, Fig. 5 shows the case of two control 
Inputs and two measurements where both the number 
of measurements and the number of active control 
bearings are less than the number of controlled 
modes. Even If there are only two active control 
bearings, the maximum effect is afforded by 
selecting the positions of measurements so as to 
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Figure 1. - Modal control system with velocity feedback. 
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Figure 4. - Eigenvalue assignment and unbalance response in the 
case of three inputs and two outputs in which J is based on the 
third mode. 
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based on the first mode* 
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